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We consider the ground state reorganization driven by an increasing nearest neighbor repulsion 
U for spinless fermions in a strongly disordered ring. When [7^0, the electrons form a glass 
with Anderson localized states. At half filling, a regular array of charges (Mott insulator) 
is pinned by the random substrate when U ^ oo. Between those two insulating limits, 
we show that there is an intermediate regime where the electron glass becomes more liquid 
before crystallizing. The liquid-like behavior of the density-density correlation function is 
accompanied by an enhancement of the persistent current. 



1 Introduction 

The interplay between the quantum interferences due to a random substrate and the correlations 
induced by charge repulsion is a central problem of mesoscopic physics. Its understanding 
seems to be a necessary step towards explaining the recently discovered metallic phase in two 
dimensions and the large value of the persistent currents observed in mesoscopic rings i. In 
one dimension (Id), assuming a Luttinger liquid in the clean limit and using renormalization 
group arguments for including elastic scattering, one finds 13 that repulsive interactions cannoJj 
delocalize spinless fermions. However, two particles in a disordered chain can be (de)localizeda 
on a length L2 larger than the one-particle localization length Li when a short range interaction 
is present. This is valid for two particle excitations of sufficient energy while a ground state built 
from two one-particle states localized far away from each other remains obviously unchanged. 
Interesting effects for the ground state require reasonably high filling factors. This led us to 
study the ground state of half filled strongly disordered rings and to find that the conclusions 
of Ref. i are no longer valid when Li becomes of the order of the separation between the 
carriers. In this case, charge repulsion gives rise to successive reorganizations of the ground 
state, accompanied bv a substantial delocalization of the particles and a large enhancement of 
the persistent currentH. The interactions Uc for which the charge density is reorganized fluctuate 
from sample to sample. Thus, the delocalization effect becomes essentially negligible, though 
visible after ensemble averaging. But it becomes very striking if one studies the orbital response 
of individual mesoscopic samples as a function of U. 

At a finite density and when [/ — > 0, the electrons form a glass with Anderson localized states 
filled up to the Fermi energy. This Fermi glass is characterized by an inhomogeneous charge 
configuration imposed by the random substrate. When [/ — > 00, at half filling for a nearest 
neighbor repulsion U, the electrons form a Mott insulator (or a pinned Wigner crystal for 
long range Coulomb repulsion and lower filling factors). It is characterized by large crystalline 
domains, possibly separated by some remaining dislocations which vanish as C/ — > 00. The 
crossover between those two different insulators yields successive spatial reorganizations of the 
ground statell which correspond to weakly avoided crossings between the ground state and the 
first excitation. At those level crossings, the system becomes more sensitive to any external 
perturbation as a twist in the boundary conditions or a flux through a ring, and this enhances 
the persistent currents. 



2 Model and method 



We consider spinless fermions on a chain of L sites with nearest neighbor (NN) interaction 

L L L 

H = ^(c|ci_i + c\_^Ci) + ^ Vim + rnui-i (1) 

i=l i=l 1=1 

and twisted boundary conditions, cq = exp(i$)cL. The operators q (c|) destroy (create) a 
particle on site i and rii = clci is the occupation operator. The on-site random energies Vi are 
drawn from a box distribution of width W. The strength of the disorder W and the interaction 
U are measured in units of the kinetic energy scale (t = 1). We study interaction between NN 
at half filling, where the ground state will be a periodic array of charges located on the even 
or odd sites of the chain when U ^ oo. The numerical results are obtained using the density 
matrix renormalization group (DMRG) algorithm Q. Ground state properties in disordered Id 
systems can be calculated with an accuracy comparable to exact diagonalization, but for much 
larger systems Q. 



3 Crossover glass-liquid-crystal for increasing repulsions 

The reorganization of the ground state induced by the NN repulsion is shown in Fig. ||, where 
the density p (expectation value of n^) is plotted as a function of U and site index i. To favor the 
inhomogeneous configuration, the disorder is taken large {W = 9) and Li ~ is of order of 

the mean spacing k^^ = 2 between the charges. For C/ssO, one can see a strongly inhomogeneous 
and sample dependent density, while for large U a periodic array of charges sets in. These two 
limits are separated by a sample dependent crossover regime. 




Figure 1: Top left: Charge configuration for a typical sample (d) for 10 particles on L = 20 sites 
at = 9. Right: Density-density correlation parameter 7 for four samples with N — 10, L — 20 and 
W — 9. Thick dots: Average over up to 166 samples. Bottom left: Average of 7 in the region around 

U = Up. 

In order to quantitatively describe this reorganization of the charge density, we calculate the 
density-density correlation function 

1 ^ 

Cir) = j;^Y.PiPi+r (2) 

i=l 



for values < r < L/2. The parameter 7 = maXr{C{r)} — min,.{C(r)} is used to distinguish 
between the electron liquid with constant density (7 = 0) and the regular crystalline array of 
charges (7 = 1). If one includes the translation r = in the definition of 7, one gets 7 7^ for 
the electron glass. So 7 measures charge crystallization from an electron liquid as well as the 
melting of the glassy state towards a more liquid ground state. 

Fig. H (right) shows the dependence of 7 on the interaction strength U for four individual 
samples. For certain impurity configurations, like in sample (a), the periodic array is obtained at 
a weak repulsive interaction, while one needs a strong interaction for other samples like (b) and 
(d). Typically, 7 assumes a minimum for a small repulsive interaction of the order of the kinetic 
energy scale t. This means that the charge distribution is closest to a liquid there and suggests a 
maximum of the mobility of the charge carriers. This is an indication for a delocalization of the 
ground state by repulsive interactions. In addition, most of the samples show small steps in the 
interval < U < 2t, caused by instabilities between different configurations of similar structure. 
The formation of the regular array of charges imposed by strong repulsive interactions occurs 
only at an interaction strength U ~ f7w with C/w oc W. It can be inferred from the step-like 
increase of 7 in individual samples that the regular array is established abruptly at an interaction 
value which depends on the disorder realization of the given sample. Therefore, the jumps of 7 
are smeared out when the average over the ensemble is calculated. 

Nevertheless, it can be clearly observed that the charge density does not cross over smoothly 
from the disordered Fermi glass to the regular Mott insulator. There is a regime at weak repul- 
sive interaction U ^ Up ^ t, where the interaction delocalizes the particles and counteracts the 
disorder before the tendency towards the Mott insulator starts to dominate at stronger interac- 
tion strength. This is similar to the behavior observed in 2d systems EpIj, where an intermediate 
phase has been found at moderate Coulomb repulsion and an abrupt transition to a regular 
Wigner crystal when the interaction strength increases. 

4 Enhanced persistent currents for intermediate repulsions 

To measure the delocalization effect associated to the change of charge configuration, we study 
the phase sensitivity of the ground state. The energy difference between periodic (<1> = 0) and 
anti-periodic (<I> = 7r) boundary conditions, AE = {—)^{E{0)—E{7r)) conveys similar information, 
in the localized regime, as other measures of the response of the ground state to an infinitesimal 
flux threading the ring: the Kohn curvature (charge stiffness) oc E ($ = 0) and the persistent 
currentij oc —£''(<!> = 0). For strictly Id systems, the sign of £'(0)— ^'(vr) simply depends on the 
parity of N, and the factor (— )^ makes AE positive. The phase sensitivity D(U) = {L/2)AE 
is shown in Fig. |^ for four samples at half filling with W = 9. Both for C/~0 and U^l, D{U) 
is very small, but sharp peaks appear at sample dependent values Uc, where D{Uc) in certain 
samples can be 4 orders of magnitude larger than for free fermions. Remarkably, the curves for 
each sample do not present any singularity at C/ = which could have allowed to locate the free 
fermion case. Peaks can be seen at different sample dependent values of U (positive or negative). 
A comparison of the peak values with the step positions in Fig. |l] shows that the sharp peaks 
are in all cases accompanied by a change of the ground state structure. 

For small repulsive interactions, the system is an Anderson insulator delocalized by U, and 
D{U) increases as a function of U, consistent with the decrease of 7 seen in Fig. |l[ This 
decrease means that one has first a progressive melting of the electron glass when one turns on 
the interaction. At the maximum of D(U) occurring at the largest value of the interaction, the 
regular array of charges is established, and thereafter it becomes more and more rigid (pinned 
by the random lattice); thus D{U) decreases as a function of U. 

The thresholds Uc are strongly sample dependent giving rise to a very wide distribution 
of D(U): the ensemble average at a given U mixes very different behaviors and provides very 




Figure 2: Left: Phase sensitivity D{U) for the four different samples of Fig. ^ in (decimal) logarithmic 
scale. Thick dots and inset: average of log(-D). Right: The energy spacing A between the ground state 
and the first excited state of sample (d) (upper line) together with the corresponding phase sensitivity 

(lower line). 



incomplete information. As shown in Fig. ^, (log !)([/)) decreases for repulsive interactions, 
except for a small interval around Uy (inset) where a local maximum is obtained. This is 
reflected in the behavior of the ensemble average (7) (Fig. 2), which increases smoothly at large 
interaction except for a maximum at U-p. 

This delocalization effect of the ground state for U ^ Up only occurs at strong disorder. For 
weak disorder {W = 2, L) we recover^ the behavior expected starting from the clean limit, 
using bosonization and renormalization group arguments 0't3: a repulsive interaction reinforces 
localization, in contrast to a (not too strong) attractive interaction which delocalizes. How- 
ever, the conclusion that a repulsive interaction favors localization is no longer valid at strong 
disorder. Similar conclusions have been drawn from a study of the conductance of one- and 
two-dimensional systems at half filling 0. 



5 Avoided level crossings 

The sharp changes in the ground state structure and the peaks observed in the phase sensitivity 
are the consequences of avoided crossings between the ground state and the first excitation, 
when U increases. While the ground state at weak interaction is well adapted to the disordered 
potential, another state with a different structure being better adapted to repulsive interactions, 
becomes the ground state at stronger interaction. In the case of large disorder, with a one-particle 
localization length Li which is of the order of the mean distance between the particles, the overlap 
matrix elements between the different noninteracting eigenstates due to the interaction are very 
small and the levels almost cross. There is only a very small interaction range where a significant 
mixing of two states is present. This is exactly where the peaks of the phase sensitivity appear. 

This scenario is confirmed by Fig. ^ (right), where the phase sensitivity and the energy level 
spacing A between the ground state and the first excited state are shown for sample (d). A 
minimum of A appears always at the same interaction value as a peak in the phase sensitivity. 
A gap of increasing size opens between the ground state and the first excitation after the last 
avoided crossing at interactions U > C/wi when the Mott insulator is established. A study of 
many other samples leads to the same conclusions. It is interesting to investigate the statistics 
of the first excitation energy. In 2djkSj'ith Coulomb repulsion, this has recently been addressed 
giving rise to intermediate statistics Ii3 at the opening of the quantum Coulomb gap. 



6 Conclusions 



We summarize our main conclusions: (i) Observables for individual samples convey a clearer 
information than the (log)-averages over the ensemble; (ii) Spinless noninteracting fermions 
in Id form an insulating Fermi glass at strong disorder, which becomes a Mott insulator when 
U C/w oc W; (iii) At an intermediate repulsion U ^ Up ^ t < C/w, the charge density is closest 
to a liquid and the mobility of the charge carriers becomes maximum, indicating the existence 
of a regime where the interaction delocalizes between two insulating limits. In small 2d clusters 
with Coulomb interaction, a rather similar behavior appliesiil. From a study of the longitudinal 
and transverse currents driven by a flux threading the system around its longitudinal axis, 
an intermediate phase was found, characterized by the suppression of the transverse currents 
and the persistence or the enhancement of the longitudinal currents. According to transport 
measurements, this new phase might correspond to a new kind of 2d rrmtal. In one dimension, 
the study of D{U) when the size L increases at a given filling factor t3 allows us to exclude 
the existence of a Id metal. We can neither characterize a change in the ground state by the 
suppression of the transverse current as in 2d. However, from a study of the density-density 
correlation functions for the ground state, we have shown that also in Id, the ordered array of 
charges is established at an interaction strength ?7w larger than the value U-p where one has a 
maximum interaction induced delocalization. We believe that the regime U ^ U-p where the 
system is nor a Fermi glass neither a Mott insulator in Id can be seen as a Id precursor of the 
probable metallic phase appearing in 2d. 
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